#include <iostream>
#include "../include/BVPs.hpp"
using namespace std;
int main()
{
    BVPs S;
    S.set_grid(64);
    //S.set_disk();
    S.solve_Dirichlet_bound();
    //S.solve_Neumann_bound();
    return 0;
}